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ABSTRACT 

This  work  was  undertaken  to  evaluate  instrumentation  and  techniques 
for  observing  the  phenomena  associated  with  the  atmospheric  re-entry,  at 
orbital  or  near  orbital  velocities,  of  NAP  systems,  with  emphasis  being 
placed  on  fuel  element  re-entry  burn-up  phenomena. 

The  uses  of  pyrometry,  radiometry,  spectrophotometry,  and  photography 
for  determining  the  temperature  of  a  re-entering  body  are  analyzed.  It  is 
shown  that  all  but  the  first  of  these  methods  can  be  used  to  obtain  temperature 
information  on  re-entering  objects  that  are  not  ablating.  This  temperature 
information  is  gained  from  measurements  on  the  shape  of  the  thermal  emission 
curve.  Methods  for  handling  the  ablating  case  are  discussed. 

The  use  of  photography  for  identifying  tracer  flares  that  are  re-entering 
the  earth's  atmosphere  at  high  velocities  was  analyzed  and  shown  to  be  feasible. 
How  photography,  in  conjunction  with  tracer  flare  re-entry  experiments,  can 
be  used  for  establishing  the  validity  of  mathematical,  re-entry,  burn-up  models 
was  also  analyzed.  The  various  associated  difficulties  and  how  they  might  be 
overcome  are  discussed  in  some  detail. 
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I.  INTRODUCTION 

This  report  summarizes  an  instrumentation  study  and  analysis  performed  by 
the  General  Technologies  Corporation  for  the  Nuclear  Power  Division  of  the 
Development  Directorate  of  the  Air  Force  Special  Weapons  Center  under  Contract 
No.  AF  29(601)-5390.  It  is  concerned  almost  exclusively  with  experiments  and 
instrumentation  that  are  related  to  the  observation  of  optical  phenomena  associated 
with  ob|ects  re-entering  the  earth's  atmosphere  at  velocities  up  to  26,000  ft/sec. 

During  the  program,  the  General  Technologies  Corporation,  in  cooperation 
with  AFSWC  personnel,  studied  the  results  of  the  over-all  NAP  systems  re-entry 
program  that  is  being  conducted  by  AFSWC  and  the  AEC*  .  These  studies  and  analy¬ 
ses  formed  the  basis  for  evaluating  the  various  possible  approaches  to  the  optical 
instrumentation  of  the  NAP  system  re-entry  tests.  The  instrumentation  studied  in¬ 
cluded  pyrometry,  spectral  photometry,  ordinary  photography  in  conjunction  with 
filters,  and  spectral  photography . 


*NAP  is  an  acronym  for  Nuclear  Auxiliary  Power. 
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II.  PROGRAM  OBJECTIVES 

The  objectives  of  the  program  were  three-fold:  (1)  to  evaluate  the  instru¬ 
mentation  being  utilized  during  NAP  system  re-entry  tests,  and,  where  possible, 
to  make  recommendations  for  changes  that  might  enhance  the  probability  of  success 
in  these  tests;  (2)  to  determine  the  reliability  of  various  instrumentation  for  the 
acquisition  of  data  on  various  re-entry  phenomena  as  specified  by  AFSWC;  and 
(3)  to  determine,  if  possible,  other  phenomena  which  should  be  measured  in  order 
to  enhance  the  utility  of  the  flight  tests'  experiments  toward  establishing  a  re-entry 
burn -up  model  for  nuclear  auxiliary  power  systems. 

How  these  objectives  were  achieved  is  shown  in  the  following  technical 
discussion  section  of  this  report. 
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III.  TECHNICAL  DISCUSSION 

This  section  of  the  report  will  first  discuss  conventional  (pyrometric, 
radiometric,  and  spectrophotometric)  and  photographic  techniques  for  deter¬ 
mining  the  temperature  of  a  re-entering  object,  and  then  proceed  with  a  dis¬ 
cussion  of  the  use  of  flaring  materials  for  determining  heat  energy  inputs  to 
re-entering  bodies . 

A.  CONVENTIONAL  TEMPERATURE  MEASURING  TECHNIQUES 

The  problem  to  be  considered  is  that  of  making  temperature  measurements 
on  objects  that  are  re-entering  the  earth's  atmosphere  at  orbital  or  near  orbital 
velocities.  The  measurements  are  to  be  made  on  instruments  that  are  located 
within  aircraft  or  in  ground  installations.  The  techniques  to  be  considered  are 
pyrometry,  radiometry,  and  spectrophotometry.  The  theoretical  bases  for  the 
operation  of  these  instruments  are  the  Wein  displacement,  Stefan-Boltzmann  radi¬ 
ation,  and  Planck  radiation  laws  of  classical  physics. 

1 .  Theoretical  Basis  of  Measurement 

In  order  to  make  temperature  measurements  on  re-entering  bodies  from 
ground  or  aircraft  installations,  it  is  necessary  to  resort  to  radiation  measurements. 

Prior  to  discussing  instrumentation  for  making  such  measurements,  it  is  appropriate 
to  first  write  down  the  theoretical  basis  for  such  measurements.  Only  three  basic 
equations  are  needed  to  provide  the  basis  for  the  analysis  of  appropriate  instrumentation. 
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They  are 


e(X,  T)  c, 

P(X'T,=IV^, 

(Planck) 

(1) 

00 

P(T)  =  JP(X,T)dX  =  et  cT4 

(Stefan-Boltzmann) 

(2) 

o 

'mT  =  0-29 

(Wein) 

(3) 

P(X,  T)  =  power  emitted  per  unit  area  by  a  body  at  absolute  temperature 
T  °K  per  unit  wavelength  interval  at  the  wavelength  X. 

P(T)  =  power  emitted  per  unit  area  from  a  body  that  is  at  the  absolute 
temperature  T°K. 

e(X,  T)  =  spectral  emittance  at  the  wavelength  X  and  temperature  T.  For 
a  black  body  (total  radiator),  e(X,  T)  =  1 . 

Cf  =  total  hemispherical  emittance  of  the  surface  area  under  consider¬ 
ation.  For  a  black  body,  ey  =  1 . 

Xm  =  the  wavelength  in  centimeters  at  which  P(X,  T)/e(X,T)a$sumes  its 
maximum  value. 

ci  >  c8,  and  c  are  constants.  For  the  c.g.s.  system  of  units: 

cj,  =  3.741  x  IO”8erg-cm*/sec,  ca  =  1 .438  cm-deg.,  and 

o  =  5.669  x  10"aerg/cm*-sec-deg4. 

Only  two  radiation  measuring  techniques  lend  themselves  to  the  situation 
under  consideration.  They  are  the  total  radiation  pyrometer  and  the  spectrophotom¬ 
eter  techniques.  The  radiometer  is  actually  a  variation  on  the  spectrophotometer. 
Consider  first  the  total  radiation  pyrometer,  which  is  based  on  the  principle  of  eq.  2. 
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2.  Total  Radiation  Pyrometer 

When  a  uniformly  heated  body  of  homogeneous  composition  and  spherical 
shape  is  viewed  at  a  great  distance,  total  hemispherical  emittance  together  with 
the  atmospheric  and  inverse -square  law  effects  will  govern  the  amount  of  power 
available  for  measurement  at  the  location  of  the  pyrometer.  It  will  have  to  be 
assumed  that  ej  can  be  employed  for  bodies  having  geometrical  shapes  that  are 
not  spherical .  If  a  lens  system  is  to  be  used  with  the  pyrometer,  then  the  energy 
under  consideration  is  that  that  is  incident  per  unit  area  per  unit  time  at  the  sur¬ 
face  of  the  lens  system.  In  order  to  write  down  the  mathematical  expression  for 
this  power,  assume  u(X)  is  the  linear  absorption  coefficient  of  the  air  intervening 
between  the  source  and  point  of  measurement.  Further,  assume  that  there  is  no 
particulate  matter  in  the  light  path  so  that  build-up  factors  caused  by  multiple 
scattering  need  not  enter  the  calculation,  and  that  diffraction  is  so  small  that  its 
effects  can  be  neglected.  Then  the  power  per  unit  area  at  a  large  distance,  R, 
from  the  source  is 


A(R)  7  c,-(\,T)e'“(X)R 

P(R'T)  -  I  dX' 


(4) 


where  A(R)  is  the  total  projected  area  of  the  source  when  viewed  at  the  distance 
R.  The  R  dependence  of  A  reflects  the  possibility  that  the  size  of  the  object  may 
change  with  distance  along  the  re-entering  trajectory.  The  factor  (2  hR8)-1  comes 
into  play  from  the  inverse  square  law  and  the  fact  that  the  projected  area  is  radi¬ 
ating  into  a  solid  angle  of  2~sterodions.  Now,  to  examine  the  integrand  of 
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equation  4  in  detail . 

If  it  is  assumed  that,  at  the  temperature  T,  e(\,  T)  is  constant  and  u(\) 
does  not  vary  significantly  with  they  can  be  factored  out  of  the  integrand  to 
give,  in  accordance  with  equation  2, 


A(R)  ^  -nR 

p(R/ T)  «  e  ?T\  or  T  = 


'  2  ~  RaP(R)  “jV*  uR/4. 

aTrT^-'J  e 


(5) 


where  u  is  the  effective  linear  attenuation  coefficient  of  the  atmosphere  for  the 
light  emitted  by  the  body  in  question.  This  is  the  equation  that  one  would  have 
to  work  with  if  the  total  radiation  pyrometer  were  to  be  used  to  try  to  measure  the 
temperature  of  the  re-entering  body.  The  assumptions  necessary  to  its  derivation 
may  be  summarized  as  follows: 

1  .  The  range  R  is  accurately  known . 

2.  The  body's  surface  area  A  is  at  a  uniform  temperature. 

3.  There  is  no  significant  quantity  of  particulate  matter  in  the  atmos¬ 
phere  to  introduce  multiple  scattering  of  the  light  emitted  by  the 
body . 

4.  Refraction  effects  can  be  neglected. 

5.  The  projected  area,  A,  of  the  re-entering  body  is  known  at  the 
time  of  measurement. 

6.  The  body,  insofar  os  its  observed  emission  characteristics  are  con¬ 
cerned,  behoves  either  as  a  grey  body  or  one  of  spherical  shape  and 
homogeneous  composition. 

7.  The  emittonce,  r.t,  is  accurately  known  under  the  conditions  of 
measurement . 

8.  The  linear  absorption  coefficient  of  air  is  essentially  constant  over 
the  spectral  range  of  interest. 
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The  first  four  of  these  assumptions  are  considered  to  be  reasonably  valid. 

The  remaining  four  cast  considerable  doubt  on  the  ability  of  total  radiation  py- 
rometry  techniques  to  gather  more  than  qualitative  information.  According  to 
assumption  5,  the  projected  area  of  surface  being  observed  must  be  accurately 
known.  Since  the  body  will,  in  general,  be  asymmetrical  and  will  suffer  con¬ 
tinuous  erosion  during  re-entry,  this  factor  alone  constitutes  a  serious  problem  for 
nonspherical  bodies. 

The  assumption  that  the  body  behaves  as  a  grey  body  may  be  a  serious  one. 
Whether  or  not  this  assumption  is  valid  can  be  reasonably  well  determined  by  labo¬ 
ratory  measurements.  Although  ^  can  readily  be  measured  in  the  laboratory  up 
to  temperatures  of  2300°K,  and  a  considerable  quantity  of  data  is  available  for 
temperatures  up  to  this,  it  is  doubtful  that  it  wouid  be  available  for  any  particular 
materials  selected  at  random.  For  materials  whose  temperatures  will  exceed  2300°K 
during  re-entry,  W'H  certainly  not  be  known  except  in  rare  cases.  Thus,  assump¬ 
tion  7,  as  well  as  assumption  6,  requires  one  to  resort  to  laboratory  measurements. 

Assumption  8,  that  u(*.)  =  constant,  is  well  known  to  be  a  poor  one  over  the 
spectral  range  that  will  be  encountered.  However,  if  assumption  6  can  be  established, 
then  the  measured  rodiation  reaching  the  pyrometer  could  be  modified  for  the  absorp¬ 
tion  by  replacing  u  in  equation  5  with  the  relation 


ao 


One  would  start  with  an  assumed  value  of  |i=  Mq,  and  solve  eq.  5  for  T  =  Te.  This 
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value  of  T  would,  in  turn,  be  substituted  into  eq.  6,  using  the  known  values  of 
U(X)  for  air  to  obtain  a  better  value  of  u=  ^  .  Substituting  ^  into  5  again  gives 
a  better  value  of  T  =  T: .  This  process  can  be  reiterated  os  many  times  as  is 
necessary,  depending  upon  the  accuracy  of  the  choice  of  uQ,  until  an  accurate 
value  of  T  is  achieved. 

In  summing  up,  it  can  be  said  that,  for  a  spherically  shaped  body  whose  total 
hemispherical  emittance  is  known,  or  a  grey  body,  an  accurate  temperature  measure¬ 
ment  can  be  made  under  idealized  conditions  with  total  radiation  pyrometers  provided 
that  spectral  radiation  from  any  associated  heated  gases  or  vapors  is  negligible. 

To  see  how  the  accuracy  of  T  depends  upon  the  accuracy  to  which  et,  A(R),  R, 
and  P(R,T)  are  known,  consider  the  following  error  analysis. 

If  F  *  F(x, . . . ,  xn),  then,  according  to  the  definition  of  total  derivative, 


dF=  l^-dxi  dx9  .. .  +  dxn- 
3xx  3xg  ®  xn 

If  the  second  form  of  equation  5  is  grouped  according  to  variables,  and  these  are 

treated  os  being  independent  variables  with  xx  =  R,  x.  *  A,  x,  =  P,  x4  ■  et,  and 

xa  *  u,  then,  from  eqs.  5  and  7, 

dT  _  /  2"\yV  P  V4  [  4  ±1 

**  '  (  -J  \A7.)  T  6  —  J  2R  4 


3F 


(7) 


3T  «  T  .  9T..-I-  3T  =  -I.iT  =  RT 
sP  TV '  3A  A  ’  3ct  3u  4 
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Substituting  these  partial  derivatives  into  eq.  7,  dividing  by  T,  and  replacing  dx; 
with  Axj  gives 


AT  1  u 

T“"[2T+T 

This  equation  clearly  illustrates  the  important  point  that  errors  introduced  in  T  by 
errors  in  Liwill  be  dependent  upon  the  range  R,  and  that  errors  in  T  due  to  errors  inR 
will  be  influenced  by  the  value  of  u.  Thus,  the  larger  the  value  of  R  or  u,  the  greater 
the  error  in  T  for  a  given  fractional  error  in  uand  R  respectively.  To  see  this,  suppose 

A  i  I 

—  =  0.01  and  u*  0.04.  Then  for  Ax:  =  0  for  x;  ¥  u, 
u  1  ' 


0.01  for  R  =  100 


(.0.1  forR  =  1000 


In  a  similar  manner,  for  Axj  =  0  for  x;  fR,  u=  0.04,  'R/R  =  0.001, 


AR  - 


AA 

A 


AP 
4  P 


A  et 


(8) 


T~ 

'0.0005  +  .001 

=  0.0015  for  R  = 

100 

AT 

T  = 

(7**7] 

AR  =  < 

.0.0005  +  .01 

=  0.0105  forR  = 

1000. 

(9b) 


Thus,  it  is  apparent  that  the  error  in  T  depends  on  the  magnitude  of  both  |jand  R  as  well 
as  the  uncertainty  in  the  values  of  these  two  parameters.  These  latter  equations  show 
that  for  large  values  of  either  uor  R,  the  temperature  information  becomes  qualitative 
in  nature. 

The  conclusion  to  be  drawn  from  the  foregoing  analysis  is  that  the  total 
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radiation  pyrometer  technique  has  doubtful  utility  even  under  the  best  set  of  con¬ 
ditions  that  could  be  reasonably  assumed  in  the  intended  application. 

3.  Spectral  Photometry  Techniques 

The  application  of  spectrophotometry  techniques  should  be  considered  under 
two  separate  conditions.  The  first  of  these  is  where  the  re-entering  body  exhibits 
the  grey  body  spectrum,  that  is  T)  is  a  constant  throughout  the  wavelength  inter¬ 
val  in  which  observations  are  necessary  for  determining  the  temperature.  The  second 
case  involves  the  non-black  body  spectral  distribution  ond,  where  a  knowledge  of 
the  spectral  emittance  is  missing,  becomes  insoluble  except  when  special  attention 
is  given  to  the  emission  spectrum  curve. 

It  is  clear  from  the  foregoing  discussion  that  the  atmospheric  absorption  must 
be  taken  into  consideration  i”  using  the  spectrophotometric  technique.  However, 
since  the  principle  involved  is  identical  to  that  just  discussed,  it  will  be  left  out  of 
the  formulation  in  order  to  preserve  simplicity  of  presentation. 

3.1  Grey  body  spectrum  cose . 

This  case  can  also  be  broken  into  two  sub-cases,  where  the  emission  peak 
alone  is  considered  and  where  two  points  away  from  the  peak  ore  considered. 

3.1.1  Emission  peak  sub-case . 

In  this  case,  eq .  3  is  the  commanding  relation.  The  relative  intensity  of 
the  entire  spectrum,  or  at  least  a  portion  of  the  spectrum  that  extends  considerably 
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beyond  the  peak  emission  wavelength,  Xm  =  0.29 T-1,  on  both  sides  of  the  peak, 
is  recorded  spectrophotometry-wise.  If  this  peak  lies  in  the  vicinity  of  the  air 
absorption  curve  where  U(  X)  is  constant,  then 

T  *  0.29  X  m  for  u(X)  =  constant.  (10) 

For  the  case  where  u(  X)  is  not  constant  in  the  region  of  the  apparent  peak,  a  re¬ 
iterative  process  (as  exemplified  in  paragraph  2,  following  eq .  6)  must  be  resorted 
to.  This  is  accomplished  as  follows.  Viewing  the  absorption  and  measured  emission 
curves,  a  value  of  Xm  =  Xmo  is  chosen,  and  T  =  T0  =  Q.29/\mo  is  calculated.  This 
value  of  T0  is  used  to  manufacture  a  new  emission  curve  P( -,  Tj)  using  the  relation 

P(Ui)  =  KMo)  e‘M(X)R.  (10a) 

Here  P(X,  T0)  is  taken  from  eq.  1 ,  with  e(X,  T)  set  to  u-ity .  This  new  curve  will 
peak  at  some  value  of  X=  X^.  Using  this  value  of  Xml  and  eq.  10a,  Tj  is  com¬ 
puted.  Using  this  value  of  T  =  Tj  ,  the  second-order  approximation  of  T  =  T?,  is 
determined  from  the  curve  generated  by  the  relation 

P(X,TP)«P(\,  Tx)e'a(^R. 

The  nth  order  approximation  of  T  =  Tn  is  obtained  from  the  curve  plotted  from  the 
relation 

P(^Tn)  =  P(X,Tn.1)e"^)R,  (10b) 

where  P(  X,  Tn  _j)  is  taken  as  having  the  functional  form  of  eq.  1  with  -~(X,T)  =  1  . 

Hence  ^  =  0.29/^.  (10c) 
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The  value  of  n  required  to  obtain  a  stationii'y  value  for  the  rnonufac luted  ernis  - 
sion  peak  will  depend  upon  the  judiciousness  with  which  the  initial  value  of 
T  ~  T0  is  chosen.  Unless  m(^)  is  varying  drastically  in  the  v'cinity  of  \_.0,  one 
should  not  need  to  use  values  of  n  greater  than  2  or  3 .  It  would  perhaps  be  useful 
to  plot  the  ratio  of  P(X,  Tn)  to  the  initially  measured  curve.  When  this  intio 
approaches  a  straight  line  in  the  vicinity  of  the  peak,  one  is  assured  that  a  suffi¬ 
ciently  large  value  of  n  has  been  achieved. 

The  assumptions  employed  to  derive  equation  10c  are  much  less  capricious 
than  those  used  in  the  derivation  of  eq .  5.  For  equation  10c,  the  assumption  is 
that  e(^,T)  =  constant  in  the  vicinity  of  X  For  eq  .  10,  the  assumptions  were 
that  £(  X,  T)  and  ute  both  constant  in  the  vicinity  ol  r.m.  For  a  given  value 
of  *-m»  *he  validity  of  the  proposition  that  M(\)  is  constant  can  be  determined  from 
atmospheric  absorption  data.  The  valid! 'y  of  the  assumprion  that  £(  X,  T)  is  o  con¬ 
stant  in  the  region  of  is  something  that  will  have  to  be  determined  by  labora¬ 
tory  experiments . 


3.1.2  Two-point  or  ratio  technique  with  •  (  \ ,  T)  known  . 

When  X)  is  not  constant,  it  is  more  convenient  to  use  the  ratio  technique. 

Here  two  separate  points  on  the  emission  curve  are  chosen,  say  at  Xj  and  X9  .  Then 
using  eq .  1  and  the  exponential  atmospheric  absorption  relation,  T  is  determined 
by  the  equation. 


PU,  T) 

ca/^sT 
e  -  1 

P(X.,  T) 

^,,T)UaJ 

ecA  T.  , 

<n) 
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This  is  essentially  the  basis  for  the  operation  of  the  two-color  pyrometer.  The 
only  assumptions  required  are  that  the  body  is  either  grey  or  ^(\,T),  e(\9,T), 
are  known  and  that  and  u{\8)  are  known. 

In  view  of  the  fact  that  materials  tend  to  become  spectrally  greyer  at 
elevated  temperatures  when  the  line  spectra  are  subtracted,  this  technique  holds 
considerable  promise  for  obtaining  quantitative  temperature  information  on  a 
wide  variety  of  re-entering  materials. 

3.2  The  non-grey  bouy  cose  without  ,T). 

It  is  now  appropriate  to  see  what  can  be  done  in  case  experiments  should 
show  that  e(X,T)  is  neither  constant  nor  known.  In  other  words,  what  can  be 
determined  if  emittance  data  is  completely  lacking.  If  the  power  emission  curve 
is  recorded  as  in  the  preceding  case,  two  points  on  this  curve  can  be  employed 
to  obtain  an  expression  involving  only  one  unknown,  namely  T,  if  they  are  appro¬ 
priately  chosen.  By  viewing  the  P(X,  T)  vs.  X  curve,  there  will  be  one  or  more 
wavelength  intervals  over  which  the  emission  curve  can  be  accurately  approxi¬ 
mated  by  a  straight  line .  Over  such  a  portion  of  the  emission  curve,  the  wave¬ 
length  derivative  of  P(X,  T)  will  be  constant.  This  suggests  the  use  of  o  Taylor 
expansion  of  r(X,T)  in  this  interval,  which  can,  in  turn,  be  substituted  into  eq.1 
to  obtain  the  temperature.  In  this  case,  since  the  temperature  is  considered  to  be 
a  constant  for  o  given  spectral  curve,  the  temperature  dependence  of  e  may  be 
dropped  from  the  notation .  To  this  end,  consider  the  following. 
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Let  - — =  P'  =  constant  through  the  interval  )D  to  1  Hence 

d\ 

d  P(X,T)  _  o  f0r  n>l  in  the  interval  (X0/  X). 
d\n 

If  the  interval  (XQ/  X)  is  small,  e(\)  can  be  expanded  about  Xq  'n  a  Taylor 


€(X) 


n  w  !  ^-Xq}8  fl 

“c'*o)  n  dX!+,“Ti  dx* 


x-x„ 


4  .  .  . 

X=Xo 


For  a  given  value  of  T,  equation  I  defines  e(X)  as 


e(*)  =  P(X,T)  *(X,T), 


where 


*(X,T) 


.  V(te*/XT.  I) 


It  follows  that 

*(XQ>  =  P(X0,T)«(X0,T), 


and  that 


d* 

Jx 


_,ws  T\  d  #(X,  T) 

-^T)  -g* - 

X=X„ 


♦  ?(^o-T) 


X=X. 


dP(X,T) 

dX 


x=x„ 


In  view  of  the  fact  that  the  higher  derivatives  of  P(X,T)  are  zero,  the  highi 
rivatives  of  e(X)  can  be  written  down  immediately  os 


dn« 

dX" 


=  nP*(X0,T)^n-1)(X0,T)  4  P(X0,T)®(")(x0,T), 

x=xrt 
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where  the  notation 


d%(X,  T) 
dXm 


T)  has  been  used.  Substituting 


equations  16  and  17  into  eq.  13  gives  e(X)  in  terms  of  the  chosen  wavelengths 
XQ,  X,  the  measured  emission  power  at  XQ,  the  measured  slope  of  the  curve  at 
Xe  and  the  unknown  T.  Substituting  this  value  of  c(ty  into  equation  1  gives  the 
desired  expression  in  T.  That  is 

K  [(X  -  Xp)  P'(y0,  T)  +  P(xa,  T)] 

X^e^AT-  i) 


P(X,  T)  =  Cl 


where 


.  ?  £±£n  fil  _j _ 

ko  "«  |_<An  HKD  J  * 


(18a) 


(18b) 


The  only  unknown  in  eq.  18  is  the  temperature  T,  hence,  the  emittance  is  eliminated. 
Thus,  even  when  the  spectral  emittance  of  the  re-entering  body  is  unknown,  its  tem¬ 
perature  can  be  computed  from  the  spectral  emission  curve  by  o  judicious  choice  of 
Xa  and  X  where  the  quantity  (X  -  XQ)  is  chosen  sufficiently  small  to  assure  the  rapid 
convergence  of  the  K-series  of  eq.  18b.  Further  analysis  shows  that  the  series  con¬ 
verges  more  rapidly  for  shorter  wavelength  regions  of  the  spectral  curve,  that  is,  on 
the  short  wavelength  side  of  the  emission  peak. 

It  is  evident  that  one  of  these  three  techniques  can  be  used  to  determine  the 
temperature  providing  that  the  spectral  emission  curve  can  be  obtained. 


3.3  Unscrambling  the  spectral  emission  curve. 

The  object  is  to  devise  a  method  whereby  the  spectral  emission  produced  by 
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the  gases  and  vapors  associated  with  the  re-entering  body  can  be  separated  from 
the  observed  spectral  emission  curve.  For  the  sake  of  preserving  clarity,  the 
emission  from  the  re-entering  body  proper  will  henceforth  be  referred  to  as  the 
thermal  emission  curve.  The  observed  spectrum  (curve)  from  a  re-entering  body 
will  contain  three  components,  namely 

1 .  The  thermal  emission  curve. 

2.  The  emission  from  the  heated  air. 

3.  The  emission  from  the  vaporized  materials. 

The  thermal  emission  curve  was  discussed  briefly  above  and  is  described 
by  equation  1  . 

Up  to  temperatures  where  vaporization  and  ablotion  of  the  re-entering 
object  become  important,  one  need  only  consider  the  optical  emission  of  the 
air  in  order  to  unscramble  the  thermal  emission  curve  from  the  observed  spectrum. 
The  theoretical  basis  for  treating  both  the  heated  air  and  heated  vapors  emanating 
from  the  surface  of  the  re-entering  object  is  the  same,  so  they  will  be  considered 
simultaneously.  The  reader  will  bear  in  mind  at  all  times  that  the  nonablative 
(or  evaporative)  case  is  much  the  simpler  of  the  two  and  is  believed  to  be  well 
within  the  present  state  of  the  art  from  both  a  theoretical  and  experimental  point 
of  view. 

The  principal  mechanisms  which  produce  light  in  the  heated  material 
associated  with  the  re-entering  body  (i  .e . ,  gases  and  vapors)  are : 
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1  .  Excitation  due  to  thermal  collisions  with  the  heated  molecules 
within  the  shockwave.  Here,  within  the  shockwave  refers  to 
the  wake  as  well  as  the  volume  between  the  shock  front  and 
boundaries  of  the  re-entering  body. 

2.  Chemical  reactions  which  emit  light. 

3.  The  optically  stimulated  transitions  produced  by  nearby  sources 
of  light . 

4.  Excitation  by  impact  with  uncharged  air  molecules. 

5.  Excitation  by  collisions  with  charged  particles. 

In  considering  the  thermal  collisions,  one  con  readily  set  an  upper  bound 
to  the  quantity  of  radiation  that  it  is  possible  to  produce  from  the  kinetic  theory 
of  gases.  It  can  readily  be  shown  from  the  kinetic  theory  of  gases  that  the  frac¬ 
tion  of  molecules  having  a  velocity,  v,  greate'  than  some  specified  velocity,  vQ, 
is  given  by 


N(> * ^  .  „0./2^irr  .-m4/2kT^ 

N 


AT 


■/sA  i 


The  quantity  in  brackets  is  just  1  -  P(x),  where  P(x)  is  the  normal  probability  func¬ 
tion  and  is  tabulated  for  all  values  of  the  argument  of  interest. 


By  setting  v0=  /2hv/m’(  where  hv  is  the  energy  of  the  photon  to  be  emitted, 
eq.  19  gives  the  fraction  of  the  molecules  that  would  possess  sufficient  energy  to 
produce  an  excited  state  from  which  radiative  transitions  might  occur.  In  view  of 
the  simplicity  of  this  equation,  one  can  readily  perform  experimental  measurements 
to  determine  the  probability  per  collision  os  a  function  of  v  that  a  collision  will 
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occur  which  produces  an  excited  state  from  which  a  radiative  transition  occurs. 
In  most  practical  cases,  this  process  of  stimulated  emission  will  be  small.  There 
is  a  considerable  quantity  of  data  to  draw  from  in  evaluating  this  phenomena. 

Most  chemical  reactions  will  produce  little  light,  and  this  can  readily 
be  determined  from  the  literature  or  simple  experiments.  Thus,  this  source  of 
light  is  no  problem. 

!f  the  intensity  of  a  single  line  in  a  spectrum  is  known,  then  the  relative 

M\* 

intensities  of  other  lines,  !(v),  can  be  calculated^  '  .  This  is  illustrated  by  the 
relation 


IU)  N_(Si)_  y,4S(a,,S,) 
l(v,)  '  N(S.)  VS(a.,*) 


(20) 


where  N(  8)  is  the  number  of  molecules  in  state  B,  and  S(a,  B)  is  the  Condon- 
Shortley  line  strength^). 

In  view  of  the  fact  that  optical  transitions  involve  Einstein's  A  and  B  co¬ 
efficients,  and  of  the  fact  that  the  A  coefficient  is  derivable  from  the  B  coefficient, 
and  vice  versa,  it  is  convenient  to  define  the  quantity  5(i,  8),  as  the  probability  of 
a  spontaneous  de-excitation  transition  from  state  (or  level)  i  to  state  (or  level)  0, 
and  S(a,  i)  as  the  stimulated  excitation  transition  probability  from  level  a  to  level 
i .  It  is  understood  that  one  is  derivable  from  the  other ,  and  that  the  number  of 
sub-states  in  a  level  is  just  the  weighted  sum  of  (2J  +  I)  times  the  number  of  levels, 
where  J  is  the  total  angular  momentum  quantum  number  of  the  level .  The  optically 
stimulated  radiation  arising  from  transitions  from  the  level  0  to  the  level  a  can  then 


Numbers  in  parentheses  represent  references  at  the  end  of  the  text . 
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be  written  down  immediately  from  the  following  considerations^ . 

Assume  that  a  is  the  ground  level,  and  the  transition  from  level  6  to  level 
ais  the  one  of  interest.  For  an  exciting  intensity  l(v),  the  total  rate  of  excita¬ 
tion  from  the  ground  level  is  simply  l(\)N(a)a(  v),  where  a(v)  is  the  linear  absorp¬ 
tion  coefficient  of  the  material  containing  the  states  a.  Hence,  the  number  of 
molecules  (or  atoms)  excited  directly  to  the  level  6  is  l(v)N(a)a(v)S(a,  0).  Since 
the  mean  life  of  a  molecule  (or  atom)  in  the  excited  state  B  is  on  the  order  of  lO””7 
to  10~8  seconds,  one  may  assume  that  the  radiation  intensity  from  the  transition 
8  to  a  is  equal  to  the  rate  of  excitation  into  state  0.  The  excitation  to  these  levels 
is  simply  l(  v)N(a)a(v)  E  S(a,  i)  where  the  summation  extends  from  i  =  0  +  1  to  some 
number  n  determined  by  the  upper  energy  limit  to  which  v  is  capable  of  producing 
excitation.  Some  of  these  will  de-excite  to  level  0and  hence  to  level  a,  producing 
the  resonant  line  frequency  v0  (actually  vG  +  Av).  Since  SO*  $)  i*  the  spontaneous 
transition  probability  from  level  i  to  level  0,  the  quantity  l(v)N(a)a(v)  E  S(a,  i)  must 
be  multiplied  by  it  to  obtain  this  contribution  to  l(vQ).  Now,  all  molecules  in  level 
Bare  postulated  to  spontaneously  decay  to  level  0.  so  that  the  total  intensity  of  fre¬ 
quency  vQ  from  optical  stimulation  is 

n 

Kv0)  =  l(v)N(a)a(v)  [s(a,  B)  +  ^  S(a,  i)S(i,  8)]  .  (21) 

i=0+l 

This  can  be  normalized  to  give  the  fraction  of  excited  atoms  that  emit  the  resonant 
frequency  by  dividing  by  E  S(a,  i),  for  i  -  8  to  n.  It  is  to  be  understood  that  the 
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intensity  of  the  resonant  frequency  represents  the  intensity  of  all  optically  pro¬ 
duced  transitions  from  level  0  to  level  a.  Thus,  v0  will  exhibit  structure  when 
either  a  or  0  have  sub-levels.  For  atomic  spectra,  the  summation  over  these  sub- 
levels  can  be  treated  as  a  single  line  for  the  most  practical  cases.  It  is  perhaps 


pertinent  to  odd  the  relation  for  ®(v).  That  is. 


me  (v-v0)*+(ras/  4n)a 


where 


8”*mV  S(a,  B) 
3 he*  2J04  I 


(23) 


Here  is  the  reciprocal  of  the  mean  life  of  state  P  for  decay  to  state  a. 

The  case  of  collision  produced  excitation  can  be  treated  in  a  quite  analo¬ 
gous  manner.  It  can  readily  be  shown  that  a  similar  argument  ieads  to  the  follow¬ 
ing  relation  for  collisional  excitations. 


I(v0)r  *  N<a) 


(24) 


where  Ofis  the  density  of  the  i  th  species  of  atom  (molecules,  ions,  electrons) 
producing  the  excitation,  fr(v)  is  the  velocity  distribution  of  these  atoms  relative 
to  the  atom  (or  molecule)  with  which  they  are  colliding,  and  Qa0(v)  is  the  cross 
section  at  the  velocity  v  for  a  transition  from  level  &  to  B.  Obviously  equation 
24  must  be  summed  over  all  r  in  order  to  give  the  total  emission  due  to  collisions 
of  the  molecules  (atoms)  in  question  with  each  of  the  r  different  species  present 
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in  the  medium. 

The  above  discussion  was  intended  to  be  rather  general ;  however,  rigorous 
derivations  of  the  formulae  presented  are  derived  in  the  references  cited. 

The  geometry  of  the  volume  of  material  has  been  neglected  (optical  thick¬ 
ness),  the  nonuniformity  of  the  thermal  conditions,  the  depletion  of  states  due  to 
excitation,  the  nonequilibrium  conditions  due  to  changing  environment  during  re¬ 
entry,  etc.,  have  all  been  neglected.  However,  one  does  not  find  it  difficult  to 
imagine  that  these  can  all  be  treated  in  either  an  exact  or  an  approximate  manner 
that  will  be  sufficiently  satisfactory  to  permit  quantitative  temperature  information 
to  be  derived. 

The  literature  und  experimental  data  that  presently  exist  are  expected  to 
yield  sufficient  information  to  permit  a  determination  of  the  feasibility  of  the 
approach  in  the  nonablative  case.  A  cursory  glance  at  the  literature  shows  that 
suitable  wave  functions  and  data  are  available  for  evaluating  the  line  strengths  of 
Qs,  Ob,  QT,  O",  Np  ,  NO,  NOa,  A,  COy,  and  for  collisions  with  electrons 
for  H,  Oa,  O,  N*,  the  one  electron  (alkali)  atoms,  and  others.  It  is  worthy  of 
note  that  Kivel  and  Bailey^  were  able  to  construct  a  reasonable  model  (neglect¬ 
ing  atomic  spectra)  for  the  emissivity  of  air  as  long  ogo  as  1957.  Furthermore,  the 
composition  of  air  in  the  shock  wave  has  been  fairly  well  established  by  Trimpi 
and  Jones^,  Witliff  and  Curtis^,  and  others. 

In  conclusion,  it  is  apparent  that  any  one  of  the  three  spectrophotometric 
techniques  discussed  (emission  peak  detection,  ratio  technique,  nongrey  body  case) 
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will  successfully  determine  the  temperature  of  a  re-entering  body  provided  that  the 
line  spectrum  can  be  subtracted  from  that  part  of  the  total  spectrum  being  viewed. 

The  analysis  on  unscrambling  the  thermal  emission  curve  from  the  total  emission  curve 
indicates  that  although  the  efforts  would  be  laborous,  a  good  approximation  of  the 
thermal  emission  curve  could  be  obtained  at  least  in  the  non-ablating  re-entry  case. 
Thus,  quantitative  temperature  information  could  be  obtained  from  either  of  the 
techniques.  For  the  ablative  case,  one  would  have  to  resort  to  experiments  to  deter¬ 
mine  the  intensities  of  the  spectral  lines  that  must  be  subtracted  from  the  thermal 
emission  curve  in  order  for  these  approaches  to  work. 

B.  PHOTOGRAPHIC  TEMPERATURE  MEASURING  TECHNIQUES 

Spectral  photographic  techniques  offer  every  advantage  that  spectrophotome¬ 
try  offers  and  suffer  from  the  some  physical  restraints.  Photography  without  spec- 
tometry  is  equivalent  to  total  radiation  pyrometry. 

When  ordinary  photographs  of  a  re-entering  body  are  taken  and  the  density 
on  the  developed  negative  is  read  with  a  densitometer,  the  signal  received  is  pro¬ 
portional  to  the  logarithym  of  the  intensity  of  the  light  to  which  the  plate  was  ex¬ 
posed,  whereas  the  total  radiation  pyrometer  gives  a  signal  which  is  proportional 
to  the  intensity  of  the  incident  light.  Hence,  the  two  methods  are  essentially 
equivalent.  Thus,  this  approach  to  measuring  the  temperature  of  o  re-entering  body 
has  very  limited  utility  and  little  to  recommend  its  use. 

When  the  negative  from  a  spectral  photograph  is  read  on  a  densitometer, 
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the  result  is  a  functional  plot  whose  amplitude  is  proportional  to  the  logarithym  of 
the  intensity  at  each  wavelength  of  the  light  that  was  incident  upon  the  plate. 

From  a  knowledge  of  the  spectral  sensitivity  of  the  photographic  plate,  its  densi¬ 
tometer  trace  can  be  converted  directly  into  an  intensity  versus  wavelength  curve 
that  is  identical  to  that  that  would  be  recorded  on  a  spectrophotometer  under  the 
same  conditions.  Therefore,  the  same  methods  (grey  body  case,  emission  peak  case, 
ratio  technique,  and  non-grey  body  case)  discussed  under  the  spectrophotometric 
method  are  applicable  .  Of  course,  limitations,  such  as  atmospheric  absorption 
and  the  spectral  emission  curve  unscrambling,  are  also  applicable. 

A  major  advantage  of  the  spectral  photographic  method  for  measuring 
temperature,  when  compared  with  the  spectrophotometric  method,  is  its  inherently 
greater  speed.  Spectrophotometers  require  from  a  few  seconds  *o  a  few  minutes 
to  generate  a  spectral  curve,  whereas  a  spectral  photograph  may  be  taken  in  a  few 
milliseconds.  This  higher  speed  assures  that  the  temperature  of  the  re-entering 
body  will  not  change  appreciably  while  the  measurement  is  being  mode. 

C.  USE  OF  FLARES  FOR  DETERMINING  HEAT  INPUT 

There  are  two  types  of  flares  that  may  be  used  to  determine  the  rate  of  heat 
input  too  re-entering  body.  One  is  the  ignition  flare,  and  the  other  is  a  "burn- 
up"  or  "tracer"  type  of  flare .  Since  the  former  of  these  renders  only  the  "single 
point"  type  of  information,  the  technical  part  of  this  discussion  will  be  limited  to 
the  tracer  flare.  Since  the  primary  purpose  of  the  flare  experiments  of  interest  to 
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AFSWC,  the  AEC,  and  their  contractors  is  to  determine  the  validity  of  certain 
mathematical,  re-entry,  burn-up  models  for  nuclear  auxiliary  power  system  fuel 
elements,  the  discussion  will  be  limited  primarily  to  this  application  of  tracer 
flares. 

The  use  of  tracer  flares  to  check  the  validity  of  re-entry,  burn-up  models 
for  fuel  elements  is  predicated  upon  the  following  proposition:  Flares  having  physi¬ 
cal  characteristics  similar  to  those  of  the  fuel  elements  can  be  constructed;  the  flare 
burn-up  profile  (altitude,  velocity,  brilliancy,  etc.)  is  calculated  using  the  same 
mathematical  model  used  to  calculate  the  burn-up  of  the  fuel  elements;  and  the  flares 
can  be  subjected  to  flight  tests  that  closely  simulate  those  to  which  the  fuel  element 
is  expected  to  be  exposed.  If  the  flare  is  observed  to  behave  in  the  manner  predicted 
by  the  model,  this  raises  the  confidence  in  the  model's  predictive  capability.  The 
objective  of  the  present  study  was  to  evaluate  the  instrumentation  techniques  that 
are  appropriate  for  observing  the  re-entry  of  such  flares  and  how  such  observations 
might  be  interpreted  to  yield  useful  information  on  the  re-entry,  burn-up  model. 

To  this  end,  it  is  appropriate  to  first  develop  an  analytical  model  for  the  phenomena 
to  be  observed  and  then  see  what  ramifications  are  required  in  order  to  apply  it  to 
an  experimental  re-entry  test. 

I .  Analytical  Model  for  Flore  Observing  Instrumentation 

In  view  of  the  facts  that  photography  is  the  most  readily  applicable  method 
for  observing  re-entering  flares,  and  that  other  techniques  which  make  use  of  the 
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radiation  emitted  by  the  re-entering  flare  can  be  reduced  to  the  principles  involved 
in  photography,  the  discussion  will  be  limited  to  photographic  techniques.  Further, 
because  the  use  of  a  movie  camera  and  a  camera  with  open  shutter  will  yield  essen¬ 
tially  the  same  information,  the  camera  shutter  will  be  assumed  to  remain  open  so  as 
to  generate  a  trace  of  the  re-entering  flare  on  the  photographic  plate. 

The  camera  will  respond  to  two  sources  of  light,  the  thermal  radiation  from  the 
heated  re-entering  flare,  P(R,T),  as  represented  by  eq.  4,  and  the  light,  l(R,  ty,  that 
is  emitted  by  the  air,  and  vaporized  and  chemicol  products  from  the  flaring  materials 
(and  the  materials  in  which  it  is  capsulated  if  it  is  a  capsulated  flare).  If  flaring  mate¬ 
rials  whose  vapor  pressures  are  high  at  relatively  low  temperatures  are  employed,  then 
the  thermal  radiation,  P(R,T),  will  be  negligibly  small.  It  will  be  assumed  herein¬ 
after  that  such  is  the  case.  To  apply  the  results  to  cases  where  the  thermal  radiation 
is  not  negligible,  all  that  one  needs  to  do  is  replace  l(R,  X)  by  I  l(R,  X)  +  P(R,  X)J  in 
the  formulae  that  follow .  With  this  in  mind,  the  density  of  the  plate  produced  by 

flare  X,  containing  m  elements,  as  viewed  by  a  densitometer,  can  be  written  as 

m  n 

D(R,X)=  log[£  X  F0  (R ,  M)  Fc(\)  Ff(  \j)  ls(Xi ,  R ,  X)]  (25) 

s=t  i  -i 


where 

F  (R,  Xj)  =  fraction  of  the  light  emitted  at  the  wavelength  Xj  that  is 
transmitted  through  the  atmosphere  the  distance  R  to  the 
camera  lense. 

Fc(Xj)  ®  fraction  of  the  light  of  wavelength  Xi  incident  at  the 

camera  that  is  transmitted  through  the  lens-filter  system 
to  the  photographic  plate . 
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F^(\|)  =  film  sensitivity  at  the  wavelength  Xj,  that  is,  the  number 

of  silver  ions  activated  per  unit  intensity  of  incident  light 
at  the  wavelength  . 

ls(RAi*X)  =  intensity  of  the  light  emitted  at  the  wavelength  X;  by  the 

m  materials  in  the  flare  at  the  point  R  along  the  trajectory. 
(Actually,  it  is  the  amount  of  light  energy  emitted  at  the 
wavelength  \j,  over  a  portion  of  the  re-entry  trajectory 
that  corresponds  to  the  slit  width  of  the  densitometer  sen¬ 
sing  head).  The  subscript  s  is  used  to  denote  the  different 
materials  in  flare  X . 

The  summation  of  i  is  to  be  taken  over  all  spectral  lines  emitted  by  each  of  the  m 
re-entering  flare  materials.  A  further  simplification  of  this  relation  can  be  made 
from  the  fact  that  the  intensities  of  all  lines  from  a  single  element  in  the  flare  can 
be  written  in  terms  of  the  intensity  of  a  single  reference  line  from  that  element^). 
If  this  reference  line  is  designated  as  *.fS  for  the  sth  element  in  the  re-entering 
flare,  then  the  intensity  from  any  line,  ^|,  from  the  sth  element  in  flare  X  can  be 
written  as 


l,(R,Xi,  X)  =  Cul(R,Xr$).  (26) 

Here  Cfo  •*  a  constant  that  has  been  tabulated^)  for  every  element,  s,  of  interest. 
It  should  be  pointed  out  here  that  the  relative  values  of  Cj,  given  in  reference  1 
are  valid  only  for  a  given  element,  s.  The  tabulated  ratio  of  two  line  intensities 
from  two  different  elements,  say  C\i/CX9,  is  not  valid  for  re-entry  applications. 

If  the  further  contraction  in  notation, 


G(R,Xj)«Fa(R,Xi)Fc(Xj)  Ff(Xj), 

is  mode,  eq.  26  can  be  substituted  into  eq.  25  to  give 
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m  n 

D(R,X)  =  log[^l(R,\rs)^G(R>Xi)Ci,]  .  (28) 

*  =  i  i=i 

Now  all  of  the  quantities  under  the  second  summation  sign  in  this  equation  can  be 
calculated  from  a  knowledge  of  the  atmospheric  attenuation  at  the  distance  R,  the 
lense  system  and  filters,  and  tabulated  spectral  line  intensities.  Therefore,  for  a 
given  value  of  s  (i  .e. ,  for  a  given  element  in  the  flare),  it  is  convenient  to  set  this 
summation  equal  to  the  known  constant  C,(R).  When  this  is  substituted  into  eq.  28, 
the  result  is 

m 

D(R,  X)  -  lag  [  £  C$(R)  l(R,  \t)].  (29) 

s*i 

The  problem  is  now  reduced  to  one  of  determining  the  intensity  of  the  reference 
line  \  for  each  of  the  m  elements  contained  in  the  flaring  material  at  each  point 
along  the  trajectory  R  . 

2.  Application  to  8urn-Up  Model  Volidificotion 

It  is  at  this  point  that  the  analysis  runs  into  difficulty.  As  mentioned  above, 
the  intensity  of  a  line  is  determined  by  optically  stimulated  transitions  (eq.  21 , 22, 
23)  and  by  impact  phenomena  (eq .  24) .  Thus,  only  such  materials  as  those  listed 
on  page  21  can  be  handled  directly  from  first  principles  without  resorting  to  experi¬ 
ments  to  determine  transition  probabilities  and  collision  cross  sections. 

This  brings  up  the  important  consequence  that,  without  a  knowledge  of 
|(R ,  \rl),  nothing  can  be  ventured  as  to  haw  much  flaring  material  should  be  used 
to  produce  a  desired  flare  intensity.  Therefore,  the  intensity  of  the  light  from  a 
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flare  alone  cannot,  in  general,  be  employed  as  a  measure  of  heat  input  at  any 
particular  time  to  a  re-entering  object  without  resorting  to  some  rather  funda¬ 
mental  experimental  work. 

These  difficulties  do  not  prohibit  the  use  of  flares  if  an  adequate  burn- 
up  model  is  available.  For  example,  if  the  burn-up  is  one  simply  of  vaporization, 
then  the  rate  of  vaporization  would  be  proportional  to  the  rate  of  heat  energy 
input  to  the  flaring  object,  and  the  length  of  time  that  the  flare  would  be  visible 
would  equal  the  time  required  for  the  heat  energy  input  to  equal  the  total  heat 
of  vaporization  of  all  of  the  flaring  material .  In  addition,  the  intensity  of  the 
light  along  the  track  would  be  expected  to  be  approximately  proportional  to  the 
product  of  the  amount  of  flaring  material  deposited  per  unit  track  length  interval 
and  the  square  of  the  velocity  of  the  flare.  Usually  an  atom  will  emit  essentially 
all  the  light  that  it  is  going  to  emit  within  a  very  short  distance  (within  the  photo¬ 
graphic  resolution)  of  the  point  of  vaporization.  Since  the  amount  of  flaring  mate¬ 
rial  per  unit  track  length  is  proportional  to  the  heat  energy  input  to  the  flare,  per 
unit  track  length,  dq/vdt,  the  intensity  at  the  position  R  along  the  trajectory  would 
be  expected  to  be  approximately  proportional  to  v(R)  dq/dt.  The  validity  of  this 
last  assumption  would  depend  upon  the  amount  of  light  emitted  by  a  vaporized  flare 
atom  being  proportional  to  its  kinetic  energy  (relative  to  the  surrounding  air  mole¬ 
cules)  at  the  time  of  evaporation . 

This  analysis  shows  that,  under  the  conditions  assumed,  both  the  duration 
and  intensity  profile  of  the  flare  trajectory  could  lend  support  to  the  validity  of  a 
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burn-up  model.  The  word  support  is  underlined  because  the  correlation  between 
the  measured  and  predicted  path  length  and  intensity  profile  does  not  prove  the 
validity  of  the  burn-up  model.  Becuase  of  the  fact  that  the  mechanism  whereby 
the  flaring  material  is  caused  to  emit  light  is  not  known,  this  correlation  merely 
shows  that  the  observation  does  not  disprove  the  model.  However,  if  the  length 
of  the  trajectory  does  not  agree  with  that  predicted  by  the  model,  this  would  not 
necessarily  disprove  the  validity  of  the  model  because  the  rate  of  evaporation  at 
which  the  flare  should  first  become  visible  and  finally  become  invisible  will  be 
unknown  unless  l(R,  \.$)  can  be  calculated. 

The  up-shot  of  this  discussion  is  that  unless  measures  are  instituted  whereby 
|(R,  \Ti)  can  be  determined,  the  use  of  tracer  flares  cannot  positively  prove  or 
disprove  a  mathematical  model  of  re-entry  burn-up  .  It  is  possible,  however, 
through  the  institution  of  appropriate  experimental  tests,  to  determine  l(R,  Xrl) . 

In  fact,  if  the  re-entry  environment  can  be  accurately  specified,  one  would  be 
able  to  calculate  l(R,\rs)  from  first  principles  for  the  hydrogen-like  elements,  such 
as  alkali  metals. 

An  important  point  should  be  brought  up  here.  If  the  length  and  position 
of  the  visible  trajectory  and  its  intensity  profile  do  correlate  well  with  that  pre¬ 
dicted,  then  one  can  make  the  assumption  that  it  is  correct  and  proceed  to  calcu¬ 
late  l(R,  \t).  With  this  information,  one  can  then  devise  a  re-entry  experiment 
that  is  significantly  different  from  the  initial  experiment  but  that  is  based  upon  the 
same  mathematical,  re-entry,  burn-up  model.  If  the  data  from  this  experiment 
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correlates  closely  with  that  predicted,  one  would  not  only  be  confident  of  the 
accuracy  of  the  mathematical  model  but  would  also  have  gained  valuable  infor¬ 
mation  on  the  reactions  which  give  rise  to  the  emitted  light. 

The  confidence  in  the  predictive  capability  of  the  model  in  fuel  element 
applications  would,  of  course,  depend  upon  how  closely  the  physical  character¬ 
istics  of  the  test  flares  approach  those  of  the  fuel  elements.  (Physical  characteris¬ 
tics  here  mean  size,  shape,  density,  chemical  properties,  mechanical  properties, 
thermophysical  properties,  etc.)  If  these  two  sets  of  characteristics  differ  markedly, 
one  has  no  assurance  that  the  mechanisms  of  failure  and  burn-up  would  be  similar 
far  the  flare  and  fuel  element. 

3.  Identification  of  Multiple  Simultaneous  Flores 

During  actual  re-entry  tests,  it  is  the  general  practice  to  have  several  flares 
re-entering  simultaneously  along  closely  spaced  trajectories.  One  would  like  to  use 
photographic  means  for  determining  which  observed  trajectory  corresponds  to  which 
flare.  To  simplify  the  analysis  without  loss  of  principle,  each  flare  will  be  assumed 
to  contain  only  one  element.  Under  these  conditions,  X  may  be  set  equal  to  s  in 
eq.  29,  which  can  then  be  rewritten  as 

D(R,  s)  *  log [l(R ,  X„)  C(R)  ]  =  log  l(R,  \n)  *  B,  (30) 

where  Bj  =  log  C»(R)  is  a  constant  that  can  be  calculated  from  the  film  sensitivity, 
camera  lens-filter  system,  the  tabulated  relative  line  intensities^)  and  l(R,  X^)  once 
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s  and  R  are  specified.  It  will  be  sufficient  to  use  only  two  simultaneous,  single 
element  flares  to  illustrate  the  principle  of  identification.  Designating  these 
two  elements  by  s  =  1 ,  and  s  =  2,  eq.  30  can  be  used  to  write  two  equations, 
namely 


D(R„1)  = 

log  l(Rx,  Xrl)  +  Bj 

(31a) 

D(R.,2)  = 

log  l(R,  ,  X^)  +  Bb 

(31b) 

where  Xri  and  Xra  are  the  chosen  reference  lines  for  flare  *1  and  *2  respectively, 
and  D(RX,  1)  and  D(Ra,  2)  are  the  observed  densities  on  the  photographic  plates  (open 
shutter  case)  that  correspond  to  the  positions  Rx  and  Ra  along  the  respective  tra¬ 
jectories.  The  subscripts  on  R  are  used  to  indicate  that  the  trajectories  of  the  two 
flares  will  not  necessarily  be  identical  Now  it  is  a  fact  that  D(RX,  1)  and  D(Ra,  2) 
will  each  start  at  zero,  rise  to  some  peak  value  and  then  decrease  to  zero  provided 
the  flares  burn  up  during  the  trajectories.  Therefore,  each  of  the  traces  on  the  photo¬ 
graph  will  contain  an  infinite  set  of  values  of  densities.  Since  the  relative  magni¬ 
tudes  of  l(Rx,  X r  )  and  l(RP,  Vra)  are  unknown  (unless  they  have  been  calculated  or 
measured  under  simulated  re-entry  conditions),  no  information  can  be  gained  from 
the  traces  that  will  identify  which  of  the  two  streaks  corresponds  to  flare  *1 .  How¬ 
ever,  if  two  photographs  are  used,  one  of  which  sees  the  light  from  one  and  only 
one  of  the  flares,  while  the  other  sees  the  light  from  both,  then  positive  identification 
is  readily  mode.  It  is  assumed  here  that  time  and  positional  information  are  available 
on  both  photographic  plates. 
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This  identification  can  be  generalized  to  n-flares.  If  n  photographs  are 
made  on  n  cameras,  each  of  which  is  filtered  to  pass  the  light  from  one  and  only 
one  flare,  then  it  is  obvious  that  each  of  the  n  flares  would  be  uniquely  identi¬ 
fied.  Alternatively,  and  more  readily  realizable,  if  n  photographs  are  taken, 
each  of  which  extinguishes  the  light  from  a  different  flare,  positive  identification 
is  again  possible.  That  is,  each  flare  is  identified  by  its  absence  on  a  particular 
plate.  It  is  important  to  note  that  no  number  of  photographs,  each  of  which  con¬ 
tains  a  trace  from  every  one  of  the  flares,  is  sufficient  to  permit  the  positive  iden¬ 
tification  of  any  flare  when  more  than  one  flare  appears  simultaneously. 

In  the  case  where  each  flare  emits  a  predominantly  different  color  of  light, 
a  single  color  photograph  which  produces  true  colors  would  be  sufficient  to  uniquely 
identify  each  flare.  Similarly,  a  snapshot  spectral  photograph  could  be  employed 
to  identify  each  flare  provided  the  spectra  from  no  two  of  the  flares  are  superimposed 
on  the  photograph. 

The  pertinent  conclusions  to  be  drawn  from  the  foregoing  discussion  are 
summarized  in  the  following  section. 
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IV.  CONCLUSIONS 


TOTAL  RADIATION  PYROMETRY  TEMPERATURE  MEASUREMENTS 

For  either  a  spherically  shaped  body  whose  total  hemispherical  emittance 
is  known,  or  a  body  that  exhibits  grey-body  emission  characteristics,  measure¬ 
ments  with  a  total  radiation  pyrometer  can  be  used  to  accurately  determine  tem¬ 
perature  provided  that  the  range  to  the  body  and  the  absorption  of  the  emitted 
light  by  the  atmosphere  intervening  the  instrument  and  object  are  small  and 
accurately  known.  The  situation  in  which  line  spectra  from  the  ablated  products 
and  the  heated  gas-cap  in  front  of  the  re-entering  object  cannot  be  neglected 
must  also  be  precluded  from  consideration  far  this  case.  The  accuracy  with  which 
the  temperature  can  be  determined  is  dependent  upon  the  magnitudes  of  the  atmos¬ 
pheric  absorption  and  the  range  of  the  object  being  viewed  as  well  as  the  inaccu¬ 
racies  in  these  two  parameters.  The  analysis  has  shown  that  the  pyrometric  technique 
has  very  limited  utility,  and  other  than  being  quite  simple,  has  little  to  recommend 
its  use. 

SPECTROPHOTOMETER  TECHNIQUE  FOR  TEMPERATURE  MEASUREMENTS 

The  spectrophotometer  is  adequate  for  measuring  the  temperature  of  a  re¬ 
entering  body,  provided  that  the  spectra  from  gases  and  vapors,  associated  with  the 
re-entering  object,  can  be  subtracted  from  the  total  spectral  emission  curve  recorded 
on  such  a  device  to  yield  the  thermal  emission  curve.  (The  latter  is  described  by 
eq .  I .)  In  non-ablating  cases,  the  spectra  that  produce  the  total  emission  curve  can 
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be  unscrambled  to  yield  the  thermal  emission  curve.  This  would  be  a  tedious 
as  well  as  a  laborous  task.  The  basic  principles  involved  seem  to  offet ,  at  piesent, 
the  only  hope  for  directly  determining  the  temperature  of  a  re-entering  body  in 
a  continuous  manner  (continuous  in  t)  by  the  institution  of  optical  measurement 
techniques.  The  spectrophotometer  technique,  however,  is  not  as  good  a  vehicle 
for  implementing  those  principles  as  is  photographic  techniques.  Where  the  radia¬ 
tion  from  the  ablation  products  becomes  significant,  it  will  be  necessary  (with 
possible  exception  of  the  hydrogen-like  atom  case)  to  perform  experiments  on  these 
products  in  order  to  establish  data  from  which  their  contribution  to  the  total  spectral 
emission  curve  can  be  determined.  This  latter  statement  is  also  true  for  photographic 
techniques. 

PHOTOGRAPHIC  METHODS  FOR  MEASURING  TEMPERATURE 

The  photographic  technique  for  measuring  the  temperature  of  a  re-entering 
object  without  spectral  information  is  equivalent  to  the  pyrometric  technique, 
whereas  spectral  phototography  is  equivalent  to  the  radiometer  and  spectrophotom¬ 
eter  techniques.  Thus,  photography  without  spectral  information  is  essentially 
useless  for  the  direct  measurement  of  the  temperature  of  a  re-entering  body.  The 
use  of  a  spectral  photograph  in  conjunction  with  a  densitometer  yields  a  plot  which 
contains  information  that  is  equivalent  to  that  from  a  spectrophotometer  plot,  the 
difference  being  that  the  densitometer  gives  the  logarithym  of  the  intensity  instead 
of  the  intensity  itself.  The  principal  advantage  of  the  spectral  photography  technique 
relative  to  the  spectrophotometer  technique  is  its  inherently  greater  speed. 


34 


TDR  63-58 


Otherwise,  it  suffers  from  the  same  shortcomings  (physical  phenomena  constraints) 
as  the  spectrophotometry  technique . 

USE  Of  TRACER  FLARES  TO  VALIDIFY  BURN-UP  MODEL 

It  is  concluded  that,  without  a  detailed  knowledge  of  the  interactions  of 
the  flaring  material  with  the  re-entry  environment  which  produce  the  light  from 
tracer  flares,  the  observations  of  the  re-entry  trajectories  can  neither  prove  nor 
disprove  the  validity  of  any  particular  mathematical  model  of  re-entry  burn-up. 

This  stems  from  the  fact  that  one  cannot  calculate  the  intensity  of  the  light  at  any 
point  along  the  trajectory  without  such  knowledge  (see  eq.  29).  For  hydrogen-like 
flare  materials  (alkali  metals),  it  is  expected  that  the  interactions  giving  rise  to 
the  emitted  light  could  be  calculated  from  a  knowledge  of  the  re-entry  environment. 
The  use  of  other  flaring  materials  will  require  the  institution  of  appropriate  experi¬ 
ments  before  one  will  be  in  a  position  to  calculate  the  magnitudes  of  the  various 
interactions  that  produce  the  light  from  the  re-entering  flare.  However,  if  the 
observed  path  length,  position,  and  intensity  profile  of  the  flaring  trajectory  all 
agree  closely  with  those  predicted  by  the  mathematical  model,  they  would  lend 
considerable  credibility  to  its  correctness. 

Furthermore,  this  correlation  could,  if  it  is  not  considered  to  be  fortuitious, 
yield  valuable  information  on  the  excitation  probabilities  of  the  flaring  atoms 
in  the  re-entry  environment.  Using  this  information,  one  could  then  devise  a  re¬ 
entry  experiment  which  would  either  prove  or  disprove  the  validity  of  the  re-entry 
model.  The  confidence  in  the  predictive  capability  of  a  model,  so  established,  in 
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its  application  to  re-entering  fuel  elements  will  depend  upon  how  closely  the 
physical  characteristics  of  the  flares  approach  those  of  the  fuel  elements. 

FLARE  IDENTIFICATION 

The  analysis  shows  that  a  series  of  black  and  white  photographs,  each 
of  which  contains  an  image  of  every  one  of  two  or  more  simultaneously  re-entering 
tracer  flares,  cannot,  in  general,  be  employed  to  positively  identify  each  of  the 
flares.  This  stems  from  the  fact  that  the  intensity  of  the  light  at  any  point  along 
the  re-entry  trajectory  cannot  be  determined  beforehand  without  resorting  to  labo¬ 
ratory  experiments  to  produce  unavailable  data.  (Alkali  metals  may  be  an  excep¬ 
tion.)  However,  it  is  pointed  out  that  if  n  photographs  are  taken  of  n  re-entering 
flare  trajectories,  each  of  which  completely  extinguishes  the  light  from  a  different 
one  of  the  flares,  each  flare  can  be  uniquely  identified.  Alternatively,  each  of 
the  n  cameras  could  pass  the  light  emitted  by  one  and  only  one  of  the  flares,  with 
each  passing  the  light  from  a  different  flare. 

Further,  if  only  a  snapshot  color  photograph,  which  reproduces  the  true 
colors,  is  taken  of  the  flares,  and  if  each  of  the  flares  emits  light  of  a  predomi¬ 
nantly  different  color,  then  the  single  photograph  will  be  sufficient  to  uniquely 
identify  each  flare.  Finally,  a  spectral  photograph  will  produce  the  same  results 
as  a  color  photograph,  providing,  of  course,  that  the  spectra  from  two  or  more 
flares  do  not  overlap  on  the  photograph . 


36 


REFERENCES 


1  .  Meggers,  Corliss,  and  Scribner,  Tables  of  Spectral  Intensities,  NBS  Mono¬ 
graph  32  (29  December  1961),  U.S.  Government  Printing  Office,  Washington, 
D  .C .  Price:  $4.00. 

2.  Condon  and  Shortley,  Theory  of  Atomic  Spectra,  Cambridge  Univ.  Press, 
London. 

3.  Massey  and  Burhop,  Electronic  and  Ionic  Impact  Phenomena,  Oxford  Univ. 
Press,  London . 

4.  J.  W.  Chamberlain,  Physics  of  the  Aurora  and  Airgiow,  Academic  Press, 
N.Y.,  N.Y.  (1961). 

5.  G.  Hertzberg,  Molecular  Spectra  and  Molecular  Structure,  Vol .  I,  Spectra 
of  Diatonic  Molecules,  Van  Nostrand,  N.Y.,  N.  Y.  (1950). 

6.  B.  Kivel  and  K.  Bailey,  Table  of  Radiation  from  High  Temperature  Air,AVCO 
Research  Report  21 . 

7.  R.  L.  Trimpi  and  R.  A.  Jones,  NASA  TR-R-63  (1960). 

8.  C  .  E .  Wittliff  and  J.  T .  Curtis,  Cal .  Rept .  No.  CAL-1 1 1  (Nov .  1961). 


37 


TDR  6  3  -  SC 


No.  Cvs 


DISTRIBUTION 


1 

1 


1 

1 


1 

1 

1 

1 

1 


1 

25 

1 

2 

2 

2 


3 


HF.ADQUARTF.  US  USAF 

USAF  Dep,  The  Inspector  General  ( AF1DI),  Norton  AFP.  Calif 
USAF  Directorate  of  Nuclear  Safety  (AFINS),  Kirtland  AFB,  NM 

MAJOR  AIR  COMMANDS 

AFSC(SCT),  Andrews  AFB,  Wash  25,  DC 
AUL,  Maxwell  AFB,  Ala 

AFSC  ORGANIZATIONS 

ASD(ASRMFP),  Wr;ght-Patterson  AFB,  Ohio 
RTD(RTN-W,  Maj  Munyon),  Bolling  AFB,  Wash  25,  DC 
SSD,  AF  Unit  Post  Office ,  Los  Angeles  45,  Calif 
(SSTRE) 

(SSZMS) 

AFMslUev  CenIRRRT/,  Hol  oman  AFB,  NM 

KIRTLAND  AFB  ORGANIZATIONS 

AFSWC  (SWEH),  Kirtland  AFB,  NM 
AFWL  ,  Kirtland  AFB,  NM 
(W'LL) 

(WLR) 
l WLDW) 

( WLDN) 

(WLDC) 

OTHER  AIR  FORCE  AGENCIES 

Director,  USAF  Project  RAND,  via:  Air  Force  Liaison 
Office,  The  RAND  Corporation,  1700  M ain  Street,  Santa 
Monica,  Calif 


38 


TDR  63-58 


DISTRIBUTION  (Cont'd) 


No.  Cys 


NAVY  ACTIVITIES 

1  Commander,  Naval  Ordnance  Laboratory,  ATTN:  Dr.  Rudlin, 

White  Oak,  Silver  Spring,  Md 

1  US  Naval  Weapons  Evaluation  Facility  (N WE F)  (Code  404), 

Kirtland  AFB ,  NM 

OTHER  DOD  ACTIVITIES 

20  Hq  Defense  Documentation  Center  for  Scientific  and  Technical 

Information  (DDC),  Comeron  Stn,  Alexandria,  Va  22314 

AEC  ACTIVITIES 

4  US  Atomic  Energy  Commission,  ATTN:  DRD/Asst  Director 

for  Compact  Systems,  Mail  Station  G- 017,  Wash  25,  DC 

4  US  Atomic  Energy  Commission,  ATTN:  DRD/Asst  Director 

for  Nuclear  Safety,  Mail  Station  G- 017,  Wash  25,  DC 

6  Sandia  Corporation,  ATTN:  Mr.  V.E.  Blake,  Sandia  Base,  NM 

1  US  Atomic  Energy  Commission,  Canoga  Park  Area  Office, 

P.O.  Box  591,  Canoga  Park,  Calif 

1  University  of  California  Lawrence  Radiation  Laboratory, 

ATTN:  Mr.  Clovis  G.  Craig,  P.O.  Box  808,  Livermore,  Calif 

OTHER 

1  OTS,  Department  of  Commerce,  Wash  25,  DC 

3  Battelle  Memorial  Institute,  ATTN:  Mr.  E.L.  Foster,  505 

King  Avenue,  Columbus  1,  Ohio 

1  Institute  of  the  Aerospace  Sciences,  Inc.  ,  2  East  64th  Street, 

New  York  21,  NY 

10  General  Dynamics/ Astronautics,  ATTN:  Mr.  E.J.  Philbin, 

P.O.  Box  1128,  Mail  Zone  596  -  2,  San  Diego  12,  Calif 

6  Stanford  Research  Institute,  Menlo  Park,  Calif 

3  Vidy a  Research  and  Development,  ATTN:  Mr.  W.J.  Fleming, 

26  26  Hanover  Street,  Palo  Alto,  Calif 


39 


TDR  63-58 


DISTRIBUTION  (Cont'd) 


No.  Cvs 


3 

3 

3 

10 

3 

3 

6 

2 

1 

1 


Astropower,  ATTN:  Dr.  A  E.  Levy-Pascal,  2968  Randolph 
Avenue,  Costa  Mesa,  Calif 

Westinghouse  E lectric  Corporation,  Astronuclear  Laboratory 
ATTN:  Mr.  Alexander  L.  Feild,  Jr.,  250  Mt.  Lebanon  Blvd.  , 
Pittsburgh,  Pa 

Armour  Research  Foundation,  ATTN:  Mr.  Jack  Hedge,  10 
West  35th  Street,  Chicago  16,  111 

Space  Information  Systems  Division,  North  American  Aviation, 
Inc.,  12214  Lakewood  Blvd.  ,  Downey,  Calif,  ATTN:  Mr.  Mark 
Morris,  Dept  468 

General  Technology  Corporation,  ATTN:  Mr.  J.D,  Graves, 
2302  Willowood  Lane,  Alexandria,  Va 

AVCO  Corporation,  ATTN:  Dr.  J.G,  Lundholm,  201  Lowell 
Street,  Wilmington,  Mass 

Aerospace  Corporation,  (Library  Tech  Documents  Group, 

ATD  61-5729  Technical  Information  Center)  P.O.  Box  95085, 
Los  Angeles,  Calif 

Aerojet-General  Nucleonics,  ATTN:  Mr.  J.P.  Lehman, 
Fostoria  Way,  San  Ramon,  Calif 

Heliodyne  Corporation,  ATTN:  Mr.  Saul  Feldman,  2365 
Westwood  Blvd.  ,  Los  Angeles  64,  Calif 

Official  Record  Copy  ( Capt.  C.E.  Brunson,  III,  WLDN-2) 


40 


